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Interferometers with atomic ensembles con- 
stitute an integral part of modern preci- 
sion metrology. However, these interferome- 
ters are fundamentally restricted by the shot 
noise limit, which can only be overcome by 
creating quantum entanglement among the 
atoms. We used spin dynamics in Bose- 
Einstein condensates to create large ensem- 
bles of up to 10 4 pair-correlated atoms with 
an interferometric sensitivity —1.61+5 ? 8 dB 
beyond the shot noise limit. Our proof-of- 
principle results point the way toward a new 
generation of atom interferometers. 

Atom interferometers exploit the wave nature of mat- 
ter, providing a unique tool for modern precision metrol- 
ogy (7, 2); most prominently, they have been used to de- 
fine the second (3, 4). However, interferometers with 
uncorrected atoms are fundamentally limited by shot 
noise (5), a remnant of the atoms' particle nature. This 
limit can be overcome by creating entanglement between 
the atoms (6). Recently, this has been demonstrated for 
spin-squeezed samples (7-72). The ultimate Heisenberg 
limit, where the sensitivity scales linearly with the num- 
ber of atoms, has been predicted for highly entangled 
states, such as NOON (73) and twin Fock states (14). 
However, they have only been created with up to five 
photons (15-19) or six ions (20, 21), but not with neu- 
tral atoms. Spin dynamics has been proposed as a pos- 
sible mechanism to turn a large fraction of an atomic 



Bose-Einstein condensate into a mixture of twin Fock 
states (22, 23). 

In an interferometric measurement, an ensemble of iV tot 
particles undergoes a sequence which distributes the par- 
ticles over two independent modes +1 and —1. The pos- 
sible many-particle states are best described by the col- 
lective spin J = N tot /2 and the population imbalance 
J z = (TV+i — 7V_i)/2, where N±i are the mode pop- 
ulations in the ±1 mode. The collective spin can be vi- 
sualized on the generalized Bloch sphere, where the z- 
coordinate represents the population imbalance J z and the 
azimuthal angle represents the relative phase between the 
two modes. The spin uncertainties can be visualized as 
the extent of the state on the sphere. For instance, un- 
entangled states spin polarized along the x-direction have 
symmetric uncertainties in relative phase and number, re- 
stricted by the uncertainty relation (Fig. 1 A). A typical 
interferometer sequence results in a rotation of the in- 
put state on the Bloch sphere by an angle 9, which de- 
pends on the observable of interest. The rotation angle is 
mapped onto the population imbalance which is finally 
measured. For unentangled input states, the measure- 
ment of the rotation angle is limited by shot noise A9 = 
<j(J z )/(y/E\d(J z )/d0\) > l/y/nN^ t , where a{J z ) is the 
standard deviation of the population imbalance and n is 
the number of independent interferometric measurements. 

A measurement beyond the shot noise limit involves a 
reduction of the output uncertainty at the expense of an 
increased uncertainty of the conjugate observable, as has 
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Figure 1: Interferometric sensitivity of an unentangled 
and a twin Fock state. The blue shaded regions repre- 
sent the uncertainty of the corresponding states. (A) Geo- 
metric representation of the sensitivity of an unentangled 
input state. A rotation on the Bloch sphere is detected 
by measuring the expectation value of J z . (B) Geometric 
representation of the sensitivity of a twin Fock input state. 
Initially, the ±1 states are equally populated (J z ) =0 and 
the relative phase is undefined; the state is represented by 
a ring around the equator. After a rotation, the ± 1 states 
are still populated equally ( J z ) =0, but the tilt of the ring 
indicates fluctuations u(J z ) increasing with the rotation 
angle. 



been realized by spin squeezing (8, 9, 11, 12). An extreme 
case is provided by the twin Fock state, where the number 
of particles in the two modes is exactly equal. Because the 
relative number uncertainty vanishes, the relative phase is 
completely undetermined. The twin Fock state is hence 
represented by an arbitrarily thin line around the equator 
of the Bloch sphere (Fig. 1 B). A rotation of this state 
maps the quantity of interest not on the expectation value 
of the output ( J z ) but rather on its fluctuations ( J z ) (24). 
The resulting sensitivity overcomes the shot noise limit, 
and is only bound by the Heisenberg limit A9 oc l/iV to t- 
Such an interferometer with twin Fock input states has 
been created for very small samples, including single (75) 



and double (16, 25) pairs of photons, and a single pair of 
9 Be+ ions (20). 

We use spin dynamics in Bose-Einstein condensates for 
the creation of up to 10 4 paired neutral atoms in two dif- 
ferent spin states. To generate these nonclassical states 
of matter, we start with a Bose-Einstein condensate with 
horizontal spin orientation (Zeeman substate in p ~ 0). In 
these ensembles, collisions may produce correlated pairs 
of atoms with spins up and down (mp = ±1) (22, 23). 
These collisions are bosonically enhanced if the output 
modes are occupied. Therefore, they act as a paramet- 
ric amplifier for a finite initial population in mp = ±1 
or for pure vacuum fluctuations (26). During the para- 
metric amplification of vacuum (27), the total number of 
atoms produced in nip = ±1 and its fluctuations increase 
exponentially with time. The conjugate variable of the to- 
tal number is the sum of the two atomic phases, whose 
fluctuations are exponentially damped (28). Furthermore, 
the number difference between mp = ±1 atoms is zero 
(without fluctuations), and hence the corresponding con- 
jugate variable, the relative phase, is fully undetermined. 
The underlying physics closely resembles that of opti- 
cal parametric down-conversion in nonlinear crystals, cur- 
rently the most important technique to generate nonclas- 
sical states of light. 

The experiments are started by creating a 87 Rb con- 
densate of 2.8 x 10 4 atoms in the hyperfine state F = 2, 
mp = in an optical dipole trap (Fig. 2 A). We initiate 
the spin dynamics at a magnetic field of 1.23 G, where 
an excited spatial mode is populated and vacuum fluctu- 
ations are amplified (29). The states F = 2, mp = ±1 
are populated by spin dynamics for an optimal duration of 
15 ms (29). Subsequently, the dipole trap is switched off 
and all three spin components are recorded by absorption 
imaging (Fig. 2 B,C). 

A test of this sensitivity requires the implementation of 
an internal-state beam splitter. To drive the transition con- 
necting the F = 2, mp = ±1 states, we employ three res- 
onant microwave pulses (Fig. 3 A). On the Bloch sphere, 
the full sequence represents a rotation around the x-axis 
by an angle 9 = r fip , where r represents the duration of 
the coupling pulse and J7r is the Rabi frequency. Figure 3 
C shows the standard deviation of the normalized popula- 
tion imbalance a{J z / J) as a function of the rotation angle 
9. It illustrates that the phase to be measured in the inter- 
ferometer is mapped onto the fluctuations of the output 
instead of the expectation value. As the picture of the ro- 
tating ring on the Bloch sphere suggests, the standard de- 
viation oscillates approximately as a(J z /J) = a |sin 9\, 
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Figure 2: Preparation and analysis of the input state. (A) At a homogeneous magnetic field, a Bose-Einstein con- 
densate in a crossed-beam optical dipole trap generates atoms in the states m,F = ±1. (B) Afterwards, the trap is 
switched off and the three Zeeman sublevels are split by a strong magnetic field gradient. (C) Finally, the three clouds 
are detected by absorption imaging and the number of atoms is counted within circular masks (white line). (D) Dis- 
tribution of the relative number of atoms in nip = ±1. The measured standard deviation of the population imbalance 
a(J z ) (orange line) is well below the shot noise limit (dashed line) and mainly limited by the detection noise (dotted 
line). The experimental result for unentangled atoms (blue line) corresponds to the combination of shot noise and 
detection noise (dash-dotted line). The shaded area indicates the error of the standard deviation. Inset: Distribution 
of the total number of atoms. The total number of atoms N tot fluctuates strongly as a result of the amplification of 
vacuum fluctuations. 



where a characterizes the interferometric contrast. Simi- 
lar to optical interferometers, the contrast relies heavily on 
an identical spatial mode of the two input states. Although 
single-mode operation can be guaranteed for small sam- 
ples in tight traps, the creation of large samples demands 
more careful consideration. Generally, spin dynamics 
populates various spatial modes (30), creating multiple 
twin Fock states without a common phase relation. After 
the beam splitter, a multi-mode state yields a number dis- 
tribution which is a convolution of many single-mode dis- 
tributions - resulting in a severe reduction of the contrast. 
We have optimized the configuration for single-mode op- 
eration (29) and reach a contrast a = 0.67(0.04), close to 
the ideal value of a = l/y/2 « 0.71. 

The symmetric beam splitter obtained for 9 = ir/2 is 
of special interest. It turns the ring on the Bloch sphere 
onto the xz-plane, leading to a maximal uncertainty of the 
population imbalance, and consequently a minimal rel- 



ative phase uncertainty. It is therefore an optimal state 
to detect the phase evolution during a Ramsey sequence. 
Figure 3 B shows the distribution of the normalized popu- 
lation imbalance. For 9 = tt/2, it assumes a characteristic 
shape (31) which reflects a high probability of detecting 
most of the atoms in one of the output ports. Ideally, either 
all even or all odd numbers should appear in the distribu- 
tion of the population imbalance J z , but this parity effect 
is not accessible within our experiments due to detection 
noise and particle losses. 

Based on the beam splitter, we demonstrate that the ex- 
perimentally created state is entangled and useful for sub- 
shot-noise interferometry. The phase estimation uncer- 
tainty A9 is inferred from the state's sensitivity to small 
rotations around an arbitrary axis in the x-y plane. It 
is deduced from error propagation according to A9 = 
A J|/ (y/n | d ( J| } /d6 1 ) for n independent repetitions of 
the interferometric measurement. In our case, we only ro- 
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Figure 3: Internal-state beam splitter for nonclassical matter waves. (A) Schematic of the beam splitter sequence. (1) 
Spin dynamics initially populates the states \F = 2, rap = ±1). (2) To couple these states, the atoms in |2, —1) are 
transferred to |1, 0) by a microwave pulse. (3) Next, a pulse of variable duration r couples the states |2, +1) and |1, 0). 
(4) Finally, atoms in 1 1 , 0) are transferred to the state 1 2 , — 1 ) to enable their independent detection. (B) Distribution of 
the normalized population imbalance for two coupling pulse durations. The strongest fluctuations are obtained for a 
coupling pulse duration corresponding to a symmetric 9 = tt/2 beam splitter. The shaded area corresponds to the ideal 
result. (C) Fluctuation of the normalized population imbalance. As expected from the representation on the Bloch 
sphere, the standard deviation a(J z /J) oscillates as |sin#| /y/2 as a function of the rotation angle 9. All quantities 
were obtained for iV tot between 3, 000 and 8, 000 atoms. 



tate the state around the x-axis, which is less susceptible 
to technical noise than a Ramsey sequence due to its short 
duration. The expectation values (J?) and ((AJ 2 ) 2 ) = 

{jf) — (J%) are shown in Figs. 4 A,B for small rotation 
angles 9. The result agrees well with the ideal case includ- 
ing a finite detection noise (29). The slight deviation for 
larger angles indicates additional technical noise, such as 
magnetic field noise or microwave intensity noise. Poly- 
nomial fits (29) allow for a precise estimate of the sensi- 
tivity (Fig. 4 C) according to the error propagation. At the 
optimal point (9 = 0.015), we reach a measurement un- 
certainty A9/A9 sn = 0.83(0.1), which is -1.6li° f dB 
below the shot noise limit A0 sn = 1/ y/n(N tot ). This re- 
sult is also — 2.515'? 8 dB below the optimal classical re- 
sult achievable with our apparatus when both shot noise 
and detection noise are considered. 

The signature of useful entanglement for sub-shot- 
noise interferometry is a Fisher information larger than 
the total number of particles T > N tot (6). The Fisher 
information also provides the ultimate achievable sensi- 
tivity via the Cramer-Rao bound A9 > 1/ \JnJ-, which 
can be reached asymptotically for a large number of mea- 
surements n. Notice, however, that here, as well as in pre- 
vious experiments on squeezing, entanglement and quan- 
tum interferometry with ultracold atoms, the total num- 



ber of particles fluctuates in different realizations; thus, 
the entanglement criteria must be generalized by replac- 
ing the fixed number of particles N tot with its statistical 
average (N lot ) (32). A lower bound for the Fisher infor- 
mation is provided by T > \d(jf) /d0\ /(A J 2 ) 2 . In 
our experiment, this lower bound is maximized for an an- 
gle of 0.015 rad, resulting from the finite detection effi- 
ciency. At this angle, we obtain F / (N lot ) > 1.45t°'2g, 
where a value greater than 1 proves entanglement. This 
entanglement can be exploited for sub-shot-noise phase 
estimation in a large variety of interferometer configura- 
tions (33), including the Ramsey interferometer. 

The presented results are mainly restricted by the de- 
tection limit of 20 atoms. A detection limit of 5 atoms 
would allow for a minimum phase estimation uncertainty 
of —13.6 dB below the shot noise limit (29). Further- 
more, the region of sub-shot-noise interferometry can be 
extended up to the whole phase domain by employing 
Maximum Likelihood or Bayesian phase estimation pro- 
tocols. The combination of large samples and quantum 
enhanced sensitivity should open exciting perspectives for 
highly sensitive measurements with a new generation of 
atom interferometers. 

Bcker et al. (34) recently reported reduced atom num- 
ber fluctuations in twin-atom beams, and the group of 
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Figure 4: Phase sensitivity of the pair-correlated state. (A) Variance of the population imbalance (a(J z ) 2 } = (J 2 ) 
for a total number iVtot between 6,400 and 7,600 atoms as a function of the rotation angle 8. The orange line is 
a quadratic fit to the data (29) and the shaded area indicates the uncertainty of the fit. Orange dashed line: the 
theoretical prediction for a detection noise of 20 atoms. (B) Same for the fourth moment of the population imbalance 
((A J 2 ) 2 ) = (Jf) — {J 2 ) 2 - (C) Sensitivity of the nonclassical state. The phase estimation uncertainty A8 with respect 
to shot noise A6> sn = 1/ \/n(N lol ) (solid orange line) is obtained from the fits in (A) and (B) and compared to the 
expected sensitivity (dashed orange line). At 9 w 0.015 rad the phase sensitivity lies — 1.6ll° dB below the shot 
noise limit (black dashed line). For a range of values of 9 it is very close to the best achievable value given by the 
Cramer-Rao lower bound (dotted orange line) (29). The experimentally measured sensitivity of an unentangled state 
(solid blue line) is measured by applying a tt/2 beam splitter to a spin polarized condensate. The technical noise of the 
beam splitter causes the measured sensitivity to lie slightly above the combination of shot noise and detection noise 
(black dash-dotted line). 



M. Oberthaler has independently used spin dynamics in 
Bose-Einstein condensates to generate atomic two-mode 
entanglement detected by a homodyning technique. 
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Supporting Online Material 

1 Magnetic field dependent spin res- 
onances 

The amplification of vacuum fluctuations is only possi- 
ble for certain magnetic fields. Figure |S"T1 A shows the 
fraction of atoms transferred to the state mp = ±1 as a 
function of the applied magnetic field, after an evolution 
time of 18 ms. The resonance structure reflects the spin 
dynamics rate of various spatial modes (30). This can be 
visualized in a simplified model, in which the effective 
trap experienced by the mp = ±1 atoms is substituted 
by a cylindrical box potential (28). The corresponding 
eigenmodes resemble the TEM laser modes. Each spatial 
mode results in a spin dynamics resonance, if its eigenen- 
ergy matches the quadratic Zeeman energy of the specific 
magnetic field. We attribute the strongest resonance at 
1.23 G to a spatial mode similar to a Hermite-Gaussian 
TEM20 laser mode (28), apparent in the density profile of 
the transferred atoms. 

The magnetic field resonances are precisely described 
by means of a Bogoliubov-de Gennes theory (30). The 
presence of unstable spin modes leads to an exponential 
amplification of the population of the mp = ±1 states, 
which depends on the magnetic field (Fig. IS II B, solid 
lines). This amplification can be triggered both by a resid- 
ual number of atoms in mp = ±1 (seed atoms) and by 
pure vacuum fluctuations (26). The dashed line displays 
the calculated ratio between the populations triggered by 
these two mechanisms. Whereas the spin dynamics at low 
magnetic fields is dominantly triggered by seed atoms, the 
reduced mode overlap leads to a vanishing influence of the 
seed at high magnetic fields. In our experiments, the spin 
dynamics is initiated at a magnetic field of 1.23 G. This 
resonance provides the strongest transfer rate and a neg- 
ligible contribution of seed atoms. Thus it ensures an ef- 
ficient amplification of almost pure vacuum fluctuations. 
Note the large separation between the two most unstable 
modes at 1 .23 G, which provides single-mode operation 
as discussed in Sec. |4] 

2 Experimental sequence and imag- 
ing 

The experimental sequence starts by preparing an ultra- 
cold thermal cloud of 87 Rb in the state F = 2, mp = 2 
in an optical dipole trap. The dipole trap is generated 




Magnetic field (G) 

Figure S 1 : Magnetic field dependence of the spin dynam- 
ics. (A) Fraction of atoms transferred to the mj? = ±1 
states versus the magnetic field after an evolution time of 
18 ms. A first resonance is observed at 0.7 G and a sec- 
ond one at 1 .23 G. (B) Numerical calculation of the am- 
plification rate of the two most unstable excitation modes 
(solid lines, left axis). The positions of the maxima re- 
flect the experimentally observed resonances. The dashed 
line (right axis) shows the ratio between the population 
triggered by vacuum fluctuations P va cuum and by residual 
seed atoms P see d at 18 ms. The calculation indicates the 
almost pure quantum character of the second resonance. 



by two crossed laser beams at a wavelength of 1064 nm. 
The atoms are transferred to the state F = 1, mp = 1 
by a resonant microwave ir pulse to avoid hyperfine- 
changing collisions during the subsequent evaporation to 
Bose-Einstein condensation (BEC). A second microwave 
pulse of variable duration transfers an adjustable num- 
ber of atoms to the state F = 2, mp = 0. To pro- 
vide a well-defined starting time for the spin dynamics 
and to avoid spurious atoms, a third microwave pulse 
transfers the atoms to the state F = 1, mp = and a 
strong magnetic field gradient removes all residual atoms 
in magnetic field sensitive states. At the same time, the 
dipole trap is adjusted to its final trapping frequencies of 
2ir x (183, 151, 122) Hz. Spin dynamics is initiated by a 
final transfer to the state F = 2, mp = 0. During the evo- 
lution time, a constant magnetic field is actively stabilized 
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to 1.23 G (with shot-to-shot fluctuations below 0.5 mG). 
The chosen magnetic field allows for the controlled popu- 
lation of a single spatial mode, where vacuum fluctuations 
are amplified (Section [T}. For some of the experiments, 
the two states F = 2, mp = ±1 are coupled at the end 
of the evolution time by a series of microwave pulses as 
shown in Fig. 3 A in the main text. 
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Figure S 2: Measurement of the number of atoms. (A) 
The effective saturation intensity is determined by prepar- 
ing atoms in the mp = ±1 states using a radio fre- 
quency pulse. These ensembles are detected with dif- 
ferent imaging beam intensities. For the correct effec- 
tive saturation intensity, the number of atoms does not 
depend on the imaging beam intensity (35). Minimizing 
the standard deviation of the mean number of atoms (in- 
set) yields 3.02(6) mW/cm 2 . (B) Typical density profile 
of the atomic clouds. Disk-shaped masks (white) are cen- 
tred around the clouds of atoms in the mp = — 1 (left) 
and mp = +1 (right) states. The number of atoms within 
each mask is summed to calculate the population of each 
cloud. The size of the masks is carefully chosen to de- 
tect the complete atom clouds and to minimize the detec- 
tion noise. The effect of averaging over single pixels is 
checked to be below 1 %. 

Finally, the dipole trap is switched off to allow for self- 
similar expansion. During a time of flight of 5 ms, the 
different Zeeman states are split by a strong magnetic 
field gradient. At a final magnetic field of 3.6 G along 



the imaging direction, the atomic density is recorded by 
resonant absorption imaging on the F = 2 — s- F' = 3 
cycling transition. A pulse of circularly polarized light 
with a duration of 70 fis and an intensity of 1.6 mW/cm 2 
is imaged onto a CCD camera. The atomic column den- 
sity is inferred from these pictures by employing an effec- 
tive saturation intensity. This effective saturation intensity 
of 3.02(6) mW/cm 2 is determined according to Ref. (35) 
as shown in Fig. IS 21 A. The accuracy of the employed 
method was checked in Ref. 9. The total number of atoms 
in each cloud is calculated by a summation of all single- 
pixel atomic densities within a disk-shaped mask (Fig 
B). 



3 Time evolution of spin dynamics 

Only for short evolution times, spin dynamics leads to the 
undisturbed parametric amplification of vacuum fluctua- 
tions. However, the creation of large samples requires a 
sufficiently long evolution time. It is hence important to 
investigate the evolution time to deduce an optimal dura- 
tion. 

The exponential amplification of the total number of 
atoms in the mp = ±1 states is shown in Fig. IS 3I A. The 
mean number increases exponentially up to an evolution 
time of 15 ms. Within this time, the standard deviation 
of the total number also increases exponentially, as pre- 
dicted for the amplification of vacuum (27). For longer 
evolution time, however, four effects lead to a saturation 
of the total number. Firstly, the number of atoms in the 
nip = condensate, initially 28, 000 atoms, is depleted, 
which reduces the production of mp = ±1 atoms. Fur- 
thermore, the density of mp = ±1 atoms is not negligi- 
ble and leads to a transfer back to the mp = state. This 
density also leads to a population of the mp = ±2 states. 
Finally, hyperfine-changing collisions result in a steady 
loss of atoms, even though the lifetime is large compared 
to the typical evolution times considered. These effects 
scale with the total number of atoms in mf = ±1 and 
therefore also reduce the standard deviation for long evo- 
lution time. 

These processes also affect the fluctuation of the popu- 
lation imbalance. Figure |S~3l B shows the standard devi- 
ation of the population imbalance cr( J z ) as a function of 
the total number of atoms iV tot for three different evolution 
times. It is compared to the detection noise of 20(1) atoms 
which we have independently measured by recording 400 
pictures without spin dynamics. Note that we define the 
detection noise in terms of J z = (iV+i — AL.i)/2. Both 
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Figure S 3: Time evolution of spin dynamics. (A) Total 
number of atoms in the m p = ± 1 states and its standard 
deviation versus evolution time. The number of trans- 
ferred atoms (circles) increase exponentially (solid line) 
depending on the evolution time. After 15 ms a saturation 
sets in (dashed line). (B) Standard deviation of the popu- 
lation imbalance a(J z ) for different evolution times. Af- 
ter 9 ms only a small number of atoms is transferred to the 
mp = ±1 states, whereas after an evolution time of 15 ms 
the number of atoms is strongly increased. For both times 
the standard deviation is dominated by the detection noise 
(dotted line). After an evolution time of 30 ms, the fluc- 
tuations of the populations of m p — + 1 and mp = — 1 
increase and approach the shot noise limit (dashed line). 
The standard deviation for a given total number of atoms 
is calculated by taking into account an interval of ±1000 
atoms. The grey shaded area indicates the error of the 
standard deviation. 



for 9 ms and 15 ms, the standard deviation is comparable 
to the detection limit and shows almost no dependence on 
the total number. For longer evolution times, represented 
here by 30 ms, the fluctuations start to rise and develop 
a dependence on the total number of atoms. These mea- 
surements result in an optimal evolution time of 15 ms, 
where a substantial population of the rriF = ±1 states is 



reached and the fluctuations of the population imbalance 
are much smaller than the detection noise. 



4 Single-mode operation 



The mp = ±1 atoms may populate multiple spatial 
modes as discussed in Section Q] This results in a mag- 
netic field dependence of the standard deviation of the 
population imbalance after the beam splitter (Fig. IS 4| A). 
For an optimal interferometric sensitivity, it is necessary 
to maximize this standard deviation. 

Insight may be gained from the simplified model us- 
ing a cylindrical box potential presented in Section \T\ 
This model shows that an optimal standard deviation is 
achieved when the most unstable spin mode is nondegen- 
erate. When this is not the case, the standard deviation is 
reduced by the square root of the number of unstable de- 
generate modes contributing to the spin dynamics, since 
the ensembles in these modes are independent. This value 
is further reduced by up to I/a/2 when spurious seed 
atoms trigger the dynamics (Fig. IS 4I B). Hence, an opti- 
mal configuration allows for the amplification of a single 
nondegenerate mode, which is significantly more unstable 
than the rest, and has a low sensitivity to spurious seeding. 

In our experiment, the spacing between low-energy 
modes of the effective trap is slightly smaller than the 
typical energy for spin-changing collisions. This leads 
to an unavoidable population of various modes, which in 
turn reduces the standard deviation of the population im- 
balance. However, during the amplification, the relative 
importance of additional modes decreases exponentially 
compared to that of the dominant one, i.e. the system ap- 
proaches the ideal single-mode operation with increasing 
time. 

A full numerical multi-mode calculation employing re- 
alistic experimental parameters shows that the resonance 
belonging to the nondegenerate TEM20 mode is most 
favourable. As expected, the standard deviation grows 
with the evolution time (Fig. IS 41 A). However, this time 
cannot be arbitrarily long as discussed in section [3] The 
best compromise is achieved for an evolution time of 
15 ms. Note that although the best value of the variance 
is achieved at very low magnetic fields around 0.4 G, it 
is unfavourable to operate at those low magnetic fields, 
mainly due to spurious seeding. 
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Figure S 4: Fluctuations of the population imbalance af- 
ter the beam splitter. (A) Numerical results for the mag- 
netic field dependence of the fluctuation of the population 
imbalance a(J z )/(J) after applying a rotation angle ir/2 
for different evolution times. The fluctuations are maxi- 
mal close to the spin dynamics resonances and present a 
maximmum at 1.3 G. (B) Fluctuations of the difference 
of the number of atoms for an evolution time of 15ms for 
various numbers of spurious seed atoms. At low magnetic 
fields fluctuations are strongly reduced by seed atoms. On 
the contrary, since the second resonance at higher mag- 
netic field is mainly triggered by vacuum fluctuations, the 
seed dependence is strongly reduced. Note that contrary 
to Fig. 3 C in the main text, the optimal standard deviation 
is one. 

5 Calculation of the interferometric 
sensitivity 

In the beam-splitting process the experimentally cre- 
ated state is transformed by the collective spin rotation 
g-jflJa, (3(,y p or a fixed value of 6, the experiment pro- 
vides a sequence of independent results, N + i and N-i, 
for the number of particles in the mp = +1 and mp = 
— 1 states, respectively. For the analysis of the phase 
uncertainty, we select the data in the interval N to t — 
N + i + iV_i = 7,000 ± 600 atoms. This is chosen for 
two reasons: First, shot noise is large with respect to the 



detection noise ^/N tot ^> 2a^ n , where Ud n ~ 20 is our 
detection noise. Second, this interval is close to the max- 
imum of the distribution of the total number of particles 
(Fig. 2 D) such that the chosen interval contains a rela- 
tively large number of data points (p ~ 100), for each 
value of 6. By using the selected data, we evaluate the 
expectation value of J" as the classical average: (J") = 

ELi «TM where ivf - (N$ - N®)/2 with 
i = 1, ...,p and a is an integer number. The values of 
(J 2 ) and (AJ 2 ) 2 = (J*) - (J 2 ) 2 obtained experimen- 
tally are shown in Figs. 4 A and 4 B for 9 values of the 
rotation angle 9. The corresponding mean square fluctua- 
tions c 2 j2} and c 2 A/2 )2, shown as error bars in Figs. 4 A 
and 4 B, are calculated from the experimental data and 
hence take the finite number of measurements and the ex- 
perimental classical noise into account. 

For a general distribution of twin Fock states, pj np = 
J2 n Pn\N/2,N/2)(N/2,N/2\ (where p N are arbitrary 
coefficients with p^ > 0, ^2 n Pn = 1), an d assuming a 
Gaussian detection noise {37), we find 

(J 2 ) gdn = (J 2 2 ) id + a 2 n (SI) 

and 

(A./, 2 )^ = (AJ 2 2 ) 2 d + 4<7^(J 2 ) id + 2ai, (S 2) 

where the expectation values are calculated for the out- 
put state p out = e _ieJx pi np e +i6IJx . The ideal expectation 
values in these expressions are 

(fa. wtim ^,, (S3) 

o 

and 

(AJ 2 ) 2 d = (l^+^-i^-^) sin^ + (J 2 ) 

(S 4) 

being (iV£ t ) = J2nPn Na - Equations <[SJ} and ( |S~2l 
are shown as dashed orange lines in Figs. 4 A and 4 B. 

The phase estimation uncertainty is obtained from the 
error propagation equation 

AJ 2 

M ~ Md(J?)/de\ ■ (S5) 

In the ideal case, by using Eqs. ( IS 3i and ( IS 4b . 
the phase estimation uncertainty strongly depends 
on the phase shift. Its optimal value, A8 = 
V^/VnVMt) + 2 (^tot>, is reached at 6 = 0. As 
shown in Fig. IS 51 the phase estimation uncertainty 
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Figure S 5: Expected phase estimation uncertainty. 
Phase estimation uncertainty with respect to shot noise 
calculated for a twin Fock input state \tp)i n p = 
\N tot /2,N tot /2) with N tot = 7000, and a Gaussian 
model of detection noise. The phase estimation uncer- 
tainty is calculated with Eq. ( IS 5b (solid lines) and the 
Cramer-Rao bound Eq. (IS 6l l (dashed lines). Different 
colors refer to different values of detection noise Od n . 
<7dn = 20 corresponds to the case of our experiment (blue 
lines). The horizontal dot-dashed line is the shot-noise 
limit. 



rapidly decreases for larger values of 
worse than the shot noise. 

For the Gaussian noise model discussed above, 
Eq. ( IS 5 b is calculated by using Eqs. ( IS II ) and (IS 21) , For 
small values of 9, the distribution of the population imbal- 
ance is dominated by the detection noise: The derivative 
d(J|) g dn/d^ vanishes, while (AJ 2 ) 2 dn remains finite. 
The phase estimation uncertainty predicted by Eq. (IS 51 
thus diverges at = 0. The minimum of Eq. ( IS 5i is ob- 
tained approximately when (J 2 ) id ~ o- 2 n . A plot of the 
phase estimation uncertainty for different values of <r& a is 
shown in Fig. IS 51 

To calculate the phase estimation uncertainty from the 
experimental results, we take error bars on the measured 
average moments into account. For each 9, we assume 
that we are in the central limit and apply a paramet- 
ric bootstrap method (38): For each experimental value 
of 9, we generate a pair of random numbers, the first 
(second) normally distributed with mean (J 2 ) [(AJ 2 ) 2 ] 

These values are then fit- 

2 , 



and variance <t 2 j 



r (AJ2) 2 l 

ted according to the function f/j2\(6) 



a\ sin 9 + b\ 



[/(aj 2 ) 2 (0) = a 2 sin 9 + 62 sin 2 9 + C2], where a\ and 
bi [0,2, 62 and C2] are fitting parameters. The choice of 
the fitting functions follows from Eqs. (IS 11 1 to (IS 41 i. The 



phase uncertainty is then calculated by using Eq. dS 51 . 
i.e. (A9) 2 = f {A j2 ) 2(9)/n\df { j2 ) (9)/d9\ 2 . This pro- 
cedure is repeated several times to guarantee statistical 
significance. The average fitting functions for (J 2 ) and 
(AJ 2 ) 2 are shown as solid orange lines in Figs. 4 A 
and 4 B, respectively. The corresponding parameters are 
01 = (6.7 ± 0.6) x 10 6 , 61 = 260 ± 100, a 2 = (3.3 ± 
0.7)xl0 13 ,6 2 = (2.5±6.6)xl0 9 ,c 2 = (1.9±0.7) x 10 6 . 
The grey shaded region in these figures indicates the mean 
square fluctuation (one standard deviation) around the av- 
erage. The average phase estimation uncertainty is shown 
as the solid orange line in Fig. 4 C. The shaded region 
in this figure is the mean square fluctuation of the phase 
uncertainty around the average value. 

For arbitrary distributions of twin Fock states, the phase 
uncertainty calculated according to Eq. (IS 5i is not opti- 
mal for all the values of 9. Quite generally, the optimal 
phase estimation uncertainty obtained by measuring the 
number of particles at the output port of the interferome- 
ter is provided by the Cramer-Rao (CR) bound 



AO, 



1 



CR 



(S 6) 



where 



and becomes T(9) = 



1 



N + i,N- 



■(- 



P(N +1 ,N^\9) 



P(N + i,N-i\6) \d9 

(S 7) 

is the Fisher information and P(iV+i, iV_i \9) are the 
conditional probabilities to measure A^+i and N-i, 
given 9. Although we do not have enough experi- 
mental data to estimate these probabilities, we can 
calculate the Fisher information within a Gaussian 
detection noise approximation which, as shown in 
Fig. 2 D, is a good model of the main source of noise 
in our apparatus. We take, as conditional probabilities, 
the convolution of the ideal i-\d(JV+i, JV_i|0) = 
(JV+i, JV_i|e- <tf J -p tap e+ <tf J -|JV +1 , JV_i) 
(39) with Gaussian distributions of 
width cr dn : P gda (N+i,N-i\6) cx 

Pid(^V + i, N_i\9). This approximation is consistent with 
Eqs. (|ST) and (ET2l . 

In the ideal case (a^ n = 0), the CR bound is 



A0 C R = 



V2 



v/(7V 2 ot ) + 2(iV tot ) 



(S 8) 



independently of the phase shift. For arbitrary small, but 
finite, detection noise, a^n > 0, the CR bound diverges 
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at = 0: In this case, the probabilities P(N +1 ,N^ 1 \9) 
are dominated by the detection noise and the derivative 
dP(iV +1 ,iV_i|0)/d0 vanishes at 9 = 0. For a dn = 20, 
the Cramer-Rao bound Eq. (IS 6l l is shown as the dotted 
orange line in Fig. 4 C. A comparison between the phase 
estimation uncertainty calculated with Eq. (IS 5i and the 
Cramer-Rao bound with Eq. (IS 6i is shown in Fig. IS 5l for 
different values of a^. It is interesting to notice that, for 
small values of 9, Eq. dS 51 l follows the CR bound, thus 
reaching the optimal uncertainty. The difference between 
Eq. dS 5i and the CR bound increases with 9 and, in par- 
ticular, is large for small detection noise. This illustrates 
the possibility to further increase the sensitivity of our ap- 
paratus by decreasing the detection noise and considering 
a more accurate phase estimation protocol, for instance, 
by taking higher moments of the distribution of the pop- 
ulation imbalance into account. Note that the main con- 
tribution to the detection noise is due to the electron shot 
noise in the camera pixels. We estimate that a reduction 
of the size of the atomic cloud by a factor of 2 (i.e. using 
shorter expansion time or Feshbach resonances {40)) and 
a camera with a quantum efficiency of 50% should allow 
for an improved detection noise cr^ w 5. 
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